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Abstract 

The distribution of the spin of the nucleon among its constituents can be paramet- 
rized in the form of polarized parton distribution functions for quarks and gluons. 
Using ah available data on the polarized structure function gi{x, Q'^), we determine 
these distributions both at leading and next -to-leading order in perturbation the- 
ory. We suggest three different, equally possible scenarios for the polarized gluon 
distribution, which is found to be only loosely constrained by current experimental 
data. We examine various possibilities of measuring polarized parton distributions 
at future experiments. 



1 Introduction 



Understanding the spin structure of the nucleon is one of the most important questions in 
strong interaction physics today. There has been renewed interest in recent years, largely 
as a result of a series of high-precision polarized deep inelastic scattering experiments. So 
far, the bulk of the information comes from measurements of the spin-dependent structure 
function gi{x, Q"^). From such measurements, spin-dependent parton distributions can be 
extracted and compared with theoretical models, or used to make predictions for future 
experiments. Several theoretical issues have received particular attention: the predictions 
of various sum-rules and the measurement of a^, the magnitude of the polarized gluon 
distribution, and the behaviour of the polarized parton distributions at small x. 

In a previous study we performed a global analysis of all available deep inelastic 
polarized structure function data in the context of the QCD-improved parton model at 
leading order. We supplemented the experimental information with theoretical assump- 
tions about the flavour content and the form of the distributions at large and small x. 
In particular, we presented three different sets of parton distributions characterized by 
qualitiatively different polarized gluon distributions, and suggested further experiments 
which could discriminate between them. 

Since the work of Ref. [|T| there have been several important advances. More pre- 
cise polarized structure function data has become available. Both the the SMC and 
SLAC-E143 collaborations |^ have measured the structure function gi of the deuteron, 
supplementing their previous proton target measurements, and the SLAC-E143 collab- 
oration has also measured the structure function g2 on proton and deuteron targets 0. 
The new gi data are particularly important for determining the polarized u- and ci-quark 
distributions separately. 

A second advance has been the theoretical calculation of the next-to-leading or- 
der contributions to the polarized splitting functions APij{x, as) which determine the 
evolution of the distributions. This not only improves the precision of the phenomenol- 
ogy, but also allows for the first time a consistent factorization/renormalization scheme 
dependence to be imposed on the analysis. 

In this study we repeat the analysis of using the new polarized structure function 
data, but now at both leading and next-to-leading order in perturbation theory. Our aim, 
as before, is to derive a consistent set of polarized parton distributions. At the same time 
we can check that the dependence seen in the data is consistent with the predictions 
of perturbative QCD. We discuss in detail the constraining power of the various data sets, 
the role of the polarized gluon distribution, and the prospects for future measurements. 
A similar analysis of polarized structure functions has been reported recently in Ref. 
using the 'dynamical parton model' approach in which the distributions at small x are 
generated dynamically from valence-like distributions at a small starting Qq scale. A 
reasonable description of the data is obtained. Our approach, in contrast, is to constrain 
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the small x distributions from the data, without prejudice as to the origin of the observed 
behaviour. This is in the same spirit as the unpolarized analysis of Refs. [|^, ^. 

We begin by reviewing the basic theoretical input to the analysis. The fundamental 
quantity of interest is the structure function gi which, in analogy with the unpolarized 
structure function Fi, can (in the QCD-improved parton model at leading order) be 
expressed in terms of the probability distributions for finding quarks with spin parallel or 
antiparallel to the longitudinally polarized parent nucleon: 

Fiix,Q') = Wel[q{x,Q')+qix,Q')] (1) 

9ii^,Q') = + (2) 

where 

? = + = - . (3) 

The first moment of gi{x,Q'^) measures the expectation value of the axial-vector 
current between two identical nucleon states: 

^ 12 36 9 ^ ' 

where the conserved nonsinglet axial vector currents are known from the hyperon decay 
constants (the EUis-Jaffe sum rule PD|): 



as = F + D 

as = 3F-D. (5) 

In contrast, the singlet axial vector current ao is non-conserved due to the axial anomaly. 
It cannot be related to any other experimental quantity, although several model-dependent 
estimates of it can be found in the literature: if the total nucleon spin is carried only by 
quarks one finds ao = 1 0, and if the strange quark sea in the nucleon is unpolarized 
then ao = ag ~ 0.579 [jlOl- Both these estimates are clearly ruled out by experimental 
measurements [g, ||, |n|, which find ao ^ 0.15 - 0.3 in the range 3 GeV^ < < 12 GeV^. 

In the QCD-improved parton model gi is expressed in terms of parton distributions 
for the polarization of quarks and gluons, 

9i{x,Q') = \M{x/y,Q')+Aq{x/y,Q') 



1 



^ (e') dy . f a,(g') 



(6) 
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The variation of these polarized parton distributions with Q is determined by the 
GLAP evolution equation |]T2|: 

d f Aq 



dlnQ^ I AG 




2n Jx y [ Eg^Pcg 




x/y,Q'). (7) 



For a consistent theoretical description of Q^), the coefficient functions in (||) 
must be truncated at the same order in perturbation theory as the splitting functions in 
(0). Although the 0{as) corrections to the coefficient functions ACq{y) and ACciy) 



Tjl have been known for a long time, a next-to-leading order study was not possible 
until very recently, as the 0{as) corrections to the APij{y) were not known. These have 
now been published for the first time in Ref. 0. 

Using the next-to-leading order evolution equations , the unpolarized parton dis- 
tributions can be determined to a high level of accuracy (see for example Refs. |0, H, |T6[) 
using a wide variety of high-precision experimental data from different processes. As 
the only data available on the spin structure of the nucleon are the structure function 
measurements of Refs. 0, ^, , we are still far from obtaining a similar level of precision 
on the polarized distributions. 

The appropriate choice of scheme for gi in the QCD-improved parton model has been 
a matter of discussion over the last few years (see for example for a recent review). 



The main issue is the treatment of the gluonic contributions to the Ellis-Jaffe sum rule. 
While the nonsinglet contributions to this sum rule can be related to conserved currents 
and are unambiguously defined, there is an arbitrariness in the decomposition of the 
singlet contribution [jl4|, depending on the absorption procedure used for the axial 
anomaly. In the MS scheme, this anomaly is absorbed into the next-to-leading order pure 



singlet quark-to-quark splitting function |19]. Another possible procedure is to allow for 



an anomalous gluonic contribution to the Ellis-Jaffe sum rule, which is absent in the 



MS-scheme. The appropriate scheme transformation ||T9[ then yields conservation of the 
quark contribution to cq. In the present study we will work entirely in the MS scheme. 
Any other choice of scheme would tie the resulting distributions to gi, as higher-order 
corrections to the coefficient functions of other processes are so far only known in the MS 
scheme .0 

The choice of the MS scheme also allows us to use our distributions in conjunction 
with a broad range of unpolarized distributions, the majority of which are defined in this 



NLO analysis of polarized deep inelastic scattering data in the so-called Adler-Bardeen scheme, 
which displays different factorization properties, has recently been carried out In Ref. the 

emphasis is on the determination of the singlet axial charge oq and the first moment of the polarized 
gluon distribution. 
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scheme. For reference, the MS scheme coefficient functions in (^) are 

1 + ^' ln^ + 2 + i/ + 5(l-y)(-2C(2)-^)| 



ACciy) = T«|2(2y-l)ln(^i-^^ +2(3-%)| (8) 

In the following section we will present the results of our leading and next-to-leading 
order parton model analyses of polarized structure function data. Section 3 discusses the 
prospects of various future experiments on polarized nucleons and Section 4 contains a 
brief summary and our conclusions. 

2 Polarized parton distributions in leading and next— 
to— leading order 

We adopt a similar approach to the global analysis of unpolarized parton distributions in 
the nucleon 0] by parametrizing the polarized distributions at the starting scale in the 
form: 

xAd^{x,Ql) = ijaAdX^'il ~ xf^il + -fdX + pdX^/^) 
xAq{x,Ql) = r]gAgx'">{l-xf^{l + -fgX + pgx^/^) 
xAG{x,Ql) = r]GAGx''^{l-xf^{l + ^GX + PGx'/'), (9) 

where we take Qq = 4 GeV^. The normalization factors are 

which ensures that the first moments of the distributions, Jq dxAf(x,Ql), are given by 
Vf- 

Various experimental measurements of unpolarized lepton(= e, p, i/)-nucleon and un- 
polarized Drell-Yan cross sections yield a reasonably precise flavour decomposition of the 
light quark {u, d, s) sea. Such a decomposition is not yet possible for the polarization of 
the light quark sea, as measurements of the structure function gi are only sensitive to 
the charge weighted sum of all quark flavours, not to the individual distributions. We 
therefore assume a SU(3)-symmetric antiquark polarization Aq{x, Ql) = Au{x, Ql) = 
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A(J(x, Qo) = ^s{x,Ql). This ad-hoc assumption is only justified at the present level of 
experimental knowledge, and is furthermore immediately broken by next-to-leading order 
evolution pi] . 



The first moments of the polarized quark distributions can be determined from the 
measured values of the Ellis-Jaffe sum rule. Imposing SU(3)-symmetry at Ql, this sum 
rule reads 

rriQl) = (l - (±^^3 + ^as + laoiQl)) (11) 

with 

^3 = rju -rjd = F + D, 
0-8 = T]u+'nd = ^F - D, 
ao(Qo) = Vu + Vd + GVgiQl)- (12) 

In this approach, the first moments of the valence quark polarizations are obtained from 
the nonsinglet axial-vector current matrix elements |2^, while the first moment of the 



sea quark distribution is inferred from the measured value of Fi. For the leading-order 
(LO) distributions, we correct the normalization of rj^ and 7]^ by the 0{as) coefficient 
function in (|TT|) []T|. The first moments obtained by this procedure are listed in Table 
Note that the Ellis-Jaffe sum rule is a (Q^-independent) constant at leading order in 
perturbation theory, as the leading-order coefficient functions are only expanded up to 
0{a^) and scaling violations in fjqiQ"^) arise only from the splitting functions at next-to- 
leading order. 

The polarized gluon distribution enters gi{x,Q'^) at next-to-leading order. It is only 
very weakly constrained so far, as no experimental data are available on gluon-initiated 
processes such as direct-7 or heavy meson production. The polarized gluon distribution 
AG{x, Q^) is therefore not well-determined by a fit to the gi data alone, and so additional 
theoretical constraints have to be applied. It follows from the structure of the polarized 
splitting functions that the small-x behaviour of the gluon and sea quark distributions are 
closely correlated, which justifies the assumption = in (P). In the region x > 0.1, 
structure functions and their evolution are dominated by valence quark contributions, and 
the impact of the gluon is completely negligible. In Ref. |jl| we explored various possibilities 
for the form of AG{x, Q^) at large x: hard and soft distributions with the spin aligned with 
that of the parent hadron, and a distribution with the spin anti-aligned. All three choices 
give equally good descriptions of the structure function data, but would be relatively 
easy to discriminate if data on polarized gluon-initiated processes were available. We 
adopt the same procedure here, i.e. we consider three equally possible scenarios for the 
behaviour of AG{x,Q'^)^ which can be parametrized as follows: 

Gluon A : 7^ = 0, pc = 0, 
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Gluon B : 7^ = -1, pa = 2, 
Gluon C : 7g = 0, Pg = —3. 



(13) 



Due to the introduction of the additional parameter pc in the starting parametrizations 
(§), these distributions look slightly different from the ones presented in |jl[. 

The normalization rjc of the gluon distribution can only be determined consistently 
from the experimental data in a next-to-leading order analysis, where it still has a large 
error. At leading order, we can estimate rjc by attributing all the violation of the EUis- 
Jaffe sum rule to a large gluon polarization and vanishing sea quark polarization. In this 
way we obtain ric = 1.9, only slightly different from the value 1.971 obtained in Ref. ||l|. 
Note, however, that at leading order the apportioning of the singlet contribution to Fi 
between gluons and sea quarks is completely arbitrary [|l^. In fact we shall see below 
that a consistent NLO treatment gives a value of 7]g somewhat less than our estimated 
leading order value, and similar to the range of values found in Ref. 
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If parton distributions are interpreted in the probabilistic picture of the naive parton 
model, the magnitude of the polarized distributions cannot exceed the unpolarized distri- 
butions, in order to guarantee positive probabilities for the individual polarization states, 
i.e. 

|A/(x)| < fix), if = q,G). (14) 

This is in fact only a rigid constraint at leading order, since parton distibutions at higher 
orders are only scheme-dependent constants of renormalization and not strict probability 
densities. The fundamental constraint at arbitrary orders in perturbation theory is the 
positivity of physical cross sections for all possible helicity configurations, which does not 
necessarily imply the positivity of the distributions. 

Positivity of the polarized distributions is achieved by constraining the parameters of 
the starting distributions at Ql such that 



\^f{x,Ql)\ < fix,Q 



0J5 



{f = q,G). 



(15) 



Perturbative evolution preserves the positivity of the individual helicity distributions, 
hence (|1^ is fulfilled at any Q^. 

In our leading-order analysis, we use the unpolarized distributions from for refer- 
ence. At Qq = 4 GeV^ these are 



xu^{x, Qq) 

xd^ (x, Qq) 

x{u + d){x, Ql) 



3.221 a;°-^^^(l - x)^-'^^^{l - 0.6889a;°-2°° + 2.254a; + 1.261a;^/2) 
^^4.476^ ^ i.615x°-^^^ + 3.651x + 1.299x^/^) 



0.507 X 



0.376/ 



^0.158 (^Q_73g _ 0.981X + 1.063a;^)(-lnx)°-°^^ + 



0.00285 exp (V-4.0101nx) (1 - x) 



\6.356 



xs{x,Ql) = 0.0034(-lnx)~^-i^(l- 2.392x^/2 + 7.094a;)(l-a;) 



6.166 
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exp (^V— 6.719 Inx ; 

xG{x,Ql) = [x°-^3i(5_iio- 1.204X- 1.911x2)(-lna;)-°-^^^^ + 

0.0527 exp (V-4.5841nx) ] (1 - xf-^^^ (16) 

The reference unpolarized distributions at next-to- leading order are the A' set of , which 
are parametrized at Qo = 4 GeV^ as 

xu^{x, Ql) = 2.26x^-^^^ (1 - xf^^ (1 - 0.54a;^/2 + 4.65a;) 

xdy{x, Ql) = 0.279a;°-23^ (1 - x)^-^^ (1 + 6.80a;^/2 + 1.93x) 

xSea(x, Ql) = 0.956x-°-^^ (1 - x)^-^^ (1 - 2.55^^/^ + 11.2x) 

xG{x, Ql) = 1.94x~°-i^ (1 - x)^-33 (1 - 1.90x^/=^ + 4.07x) (17) 

Note that the choice of unpolarized distributions is not crucial for the present analysis. 
All the widely available leading and next-to-leading order distributions provide very good 
fits to the unpolarized structure function data, and the small differences between them 
are much smaller than the precision with which the polarized distributions are currently 
determined. The flavour decomposition of the unpolarized sea quark distributions is also 
unimportant for our present analysis and will therefore be disregarded. The starting sea 
quark distribution of Ref. |^ contains a very small charm quark contribution which can 
safely be ignored in the present analysis. To be consistent with the evolution of the 
unpolarized distributions, we take 

Al^T^ = 200 MeV [16], A^^q^ = 231 MeV 0. (18) 

These correspond to as(M|) = 0.123 (LO) and as(M|) = 0.112 (NLO). 

The parameters most affected by the positivity constraints are the large-x exponents 
bf. For the valence quarks, we fix = 6„(unpol.) and b^ = ^^(unpol.) + 1, motivated by 
counting rule estimates |2^. The parameters be and bg are constrained to be at least as 



large as their unpolarized counterparts in the fit, but this constraint has only minimal 
impact. We find that only the Ad^^x) distribution tends to saturate positivity, requiring 
the combination 7^ + pd to be limited in the fit. 

The data currently available on gi are not able to test the various theoretical model 



predictions for the small- a; behaviour of the polarized parton distributions p4[ , which are 
only expected to apply at much lower values of x [^. The parameters af are therefore 
only effective exponents valid over some finite interval in x. It therefore makes no sense 
to postulate positivity for x — > by constraining the a/. 

The contribution of charmed quarks to gi{x, Q"^) is negligible at present experimental 
energies and will not be considered in this analysis. We therefore adopt the evolution 



procedure of Ref. [|T6[ and fix the number of flavours in the splitting functions at n/ = 3, 



while the number of flavours in as increases at each mass threshold, 

me = 1.5 GeV, rUb = 4.5 GeV, rrit = 180 GeV, (19) 
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LO 


NLO 


Vu 


0.823 


0.918 




-0.303 


-0.339 


Vg 


1.9 




Vq 


-0.0495 


-0.060 




3.73 


3.96 


bd 


4.73 


4.96 



Table 1: Fixed parameters in the LO and NLO (MS) polarized parton distribution fits. 



and A(n/) is determined by requiring as to be continuous across each threshold. 

Rather than measuring gi{x,Q'^) directly from absolute cross section differences, it is 
the relative asymmetry 

which is determined experimentally. The structure function gi{x,Q'^) is then inferred 
using a particular parametrization of Some experimental groups assume Q^- 

scaling of Ai{x,Q'^) in their extraction of gi{x,Q'^). In order to have a consistent set of 
data, we have instead used the measured values of Aii^x^Q"^) quoted by the experiments 
and re-evaluated gii^x^Q"^) from Eq. (pUD, constructing Q^) from the parametriza- 
tions of F2 |27] and R ||28| which were used in the most recent measurements. 

Applying the constraints outlined above, we have used all available world data on 
^^''^'"'(a;, Q^) 0, 1^, |l^ to fit the polarized quark and gluon distributions with the paramet- 
ric forms of @ imposed at Qg = 4 GeV^. About 35% of these data were taken at < Qq. 
To have sensible constraints on the distributions, in particular for x < 0.02, we include 
these datapoints in the global fit. The distributions in the region 1 GeV^ < < are 
obtained by inverting the evolution matrix, which is straightforward in A^-moment space. 

A problem with using low data points in the fit is the possible contamination by 
higher-twist contributions. We have tried to estimate to magnitude of such contribu- 
tions to gi using the parametrization of from Ref. |2^ and assuming gi{x,Q'^)^^ ~ 
(5^)^^(1 + C^'^ {x) / Q"^) . The higher-twist contributions estimated in this way are 
found to be small for all data-points apart from the two lowest x bins of the SMC exper- 
iment . 

The global fit is performed using GLAP evolution algorithms in A^-moment space |3^ . 
The distributions and structure functions are then restored by a numerical inversion into 
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X space. The results of the global fit using the leading and next-to-leading order (MS) 
expressions for the splitting functions and the gi coefficient functions are listed in Table |^. 
The resulting distributions at Qq are shown in Figs. 1 (LO) and 2 (NLO). 





A (LO) 


B (LO) 


C (LO) 


A (NLO) 


B(NLO) 


C (NLO) 


au 


0.578 


0.585 


0.582 


0.512 


0.504 


0.471 




9.38 


9.31 


9.50 


11.65 


11.98 


13.14 


Pu 


-4.26 


-4.28 


-4.28 


-4.60 


-4.61 


-4.90 


ad 


0.666 


0.662 


0.660 


0.780 


0.777 


0.809 


Id 


10.46 


10.91 


11.04 


7.81 


8.18 


6.73 


Pd 


-5.10 


-5.09 


-5.06 


-3.48 


-3.61 


-1.99 


Vg 








1.71 


1.63 


1.02 


ac 


0.520 


0.524 


0.456 


0.724 


0.670 


0.425 


be 


9.45 


6.87 


8.72 


5.71 


5.34 


11.05 


h 


15.06 


15.96 


11.82 


14.40 


18.06 


16.40 




2.30 


2.42 


2.11 


4.63 


5.30 


-2.67 


Pg 


-2.00 


-2.00 


-1.95 


-4.96 


-5.25 


-3.08 


x' 


98.3 


97.7 


100.0 


89.7 


91.0 


93.4 



Table 2: Fitted parameters in the LO and NLO (MS) polarized parton distributions at 
Qo- The values are with respect to the 110 data points included in the global fit. 

The resulting parameters are not independent of each other. In particular, a„, ad and 
ac = ttq are strongly correlated. The of the valence distributions are anticorrelated 
with the corresponding 7j, refiecting the fact that is only an effective exponent for a 
finite range in x. The 7^ and pi are also anticorrelated. The distribution is very flat 
around the local minima found by the global flts, especially with respect to the gluon and 
sea quark parameters. 

The three gluon scenarios give flts of almost identical quality, reflecting the small 
impact of the gluon distribution on gi{x,Q'^) at large and medium x. The obtained 
in the NLO flts are systematically lower due to the additional degree of freedom given by 
the normalization of the polarized gluon distribution. All flts give very good descriptions 



9 



for the polarized structure functions g^'^' {x,Q'^). This is illustrated in Fig. 3, which 
shows the NLO description of the various gi measurements using Gluon A. The curves 
correspond to = 1,4, 10,50 GeV^, reflecting the spread in values of the different 
data sets. There is a systematic decrease in the values of the data points from large 
X to small X. 

The contributions of Uv{x, Q^) and d^^x, Q^) to the neutron structure function gi{x, Q^) 
are almost equal in magnitude but opposite in sign. The neutron structure function is 
therefore much more sensitive to the sea quark polarization than gi{x, Q^) and gf{x, Q^). 
It displays a clear double peak structure, as the sea quarks are dominant in a different 
x-region than the valence quarks. A precision measurement of ^'"(a;, Q"^) will therefore 
be able to provide vital new information on the shape of the sea quark polarization. 

From a consideration of the size of the errors on the various fitted parameters, it is 
apparent that the world data on gi{x,Q'^) really only constrain the polarized valence 
quark distributions and, to a lesser extent, the overall magnitude of the sea quark po- 
larization. The flavour decomposition of the polarized sea is still completely unknown. 
Only dedicated experiments, such as the production of Drell-Yan lepton pairs or the 
flavour-tagging of final-state hadrons in polarized deep inelastic scattering, will be able 
to provide further information. Most important of all, the polarized gluon distribution 
is almost completely undetermined, as its impact on the polarized structure function is 
less than the present experimental accuracy. The variation between our three gluon sets 
certainly underestimates the true uncertainty in the distributions. 



3 Prospects for future measurements 

There are several measurements which could be feasible at future experiments and which 
yield additional information on the polarized gluon distribution. In this section we present 
some representative predictions for two of these: the dependence of gi{x, Q^), and an 
asymmetry in J/ ■i/'-photoproduction. We also briefly discuss the prospects of a measure- 
ment of polarized structure functions at the HERA collider. 

The derivative of F2{x, Q"^) with respect to has been used to measure the unpolar- 
ized gluon distribution in fixed target experiments |^2|] (0.008 < a; < 0.5) and at HERA 



(2 ■ 10^"^ < a; < 3 ■ 10^^). The method is particularly powerful at small x, where the 
gluon distribution dominates the evolution, dF2{x) / dlnQ"^ ~ Pgciv) ®G{x/y). In the 
same way, we can use our three sets of distributions A, B and C to explore the sensitivity 
of the polarized structure function evolution to AG(x, Q^). Figure 4 shows the predictions 
for the asymmetry Ai{x,Q'^) as a function of in the kinematic range representative 
of current fixed-target experiments. The asymmetry is obtained from Eq. ( PO]) with gi 
calculated using our polarized distributions and Fi calculated using the NLO unpolar- 
ized MRS(A') distributions of Ref. P. The latter are extrapolated to lower values of 
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Q^, which reproduces the full backwards evolution to within a few per cent. At large x, 
there is no sensitivity to AG{x, Q"^) - the evolution is completely dominated by the quark 
contribution. At small x, on the other hand, we see some dependence on the gluon. In 
this region AP^cd/) ® AG{x/y) < 0, and so the derivative dAi/dlnQ"^ is more negative 
for the sets that have a larger gluon polarization above the x value considered. In par- 
ticular, we see that at x ~ 0.01 the proton asymmetry is almost independent for Set 
C, but decreases with increasing Q"^ for Sets A and B. Unfortunately, the sensitivity of 
the present experiments is much worse than the differences between the various sets. To 
illustrate this, we have included two data points at a; = 0.035 from the recent SMC and 
E143 measurements. For the deuteron, the quark contribution to the structure function 
is smaller, and so the dependence on AG{x, Q^) at small x is somewhat enhanced. Con- 
sidering the large errors on the present data, it seems doubtful that a measurement of the 
polarized gluon distribution from the variation of Ai{x, Q^) is feasible for values of 
where perturbative expressions can be safely applied. Future high statistics experiments 
at higher beam energies ||3^ will clearly improve this measurement. It still has to be 
kept in mind that a determination of AG{x,Q^) from the dependence of Ai{x,Q'^) 
can never reach the quality of the corresponding unpolarized measurement, as the gluonic 
contribution is not as dominant in the evolution of polarized parton distributions as it is 
in the unpolarized evolution ||25[ . 



Inelastic J/^/j production from a nucleon target directly probes the gluon distribution 
via the photon-gluon fusion subprocess 7* + (? — >■ (cc) + g, which can be described in 
the colour-singlet model of Ref . . Depending on the virtuality of the photon, one can 
distinguish two different classes of events: photoproduction (Q^ ~ 0) and leptoproduction 



{Q"^ ^ 0). The EMC and NMC experiments have obtained measurements of 



the unpolarized leading-order gluon distribution in the range x ~ 0.05 — 0.25 from this 
process. The results agree well with gluon distributions extracted from other processes, for 
example large-pr direct photon production. The corresponding cross section for polarized 
J/ip leptoproduction has been calculated in Ref. [^. Taking the photoproduction limit 
0, one obtains 



dp'^dz 



27r 



a 



r]Gx{r],M\ 



X [A{z) + h\C{z)] ■ J^{pI, z) 



(21) 



with 



A{z) 



zHl 



(p't + (1 



M 



ZS. 



7Af 



)2 + (l-.)^(Mj/^+ 
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(1 - z)s,Nf + {s,N - M]/^f) 
S'yN — 2M]\fEj 



V 



z{l — z)s^N 



z = — — 

X = gluon helicity 
h = lepton helicity 

The polarized (unpolarized) cross sections are obtained by taking the difference (sum) over 
the helicity states. The polarized cross section is then proportional to AG{ri) = G+ijf) — 
G-{rj) and depends only on C{z), while the unpolarized cross section is proportional to 
G{rj) = G+ijf) + G^irf) and depends only on A{z). Taking the cross sections integrated 
over one of the variables, one can define two physics asymmetries: 

^,,) . 1^.^ , (22) 

MA) = ^^^^ f^-!^ . (23) 



dp^ I dp^ 

With a photon beam of energy = 45 GeV and cuts on p"^ and z (p^ > p'^^^^ = 
0.25 GeV^, z < Zmax = 0.9, as proposed in |^), one gets significantly different asymme- 
tries for the three gluon distributions A, B and C. Predictions for the asymmetry A{z), 
calculated using the three leading order distributions in combination with the unpolar- 
ized leading order gluon distributions of Ref. |T^, are shown in Fig. 5. Note that small z 
corresponds to large rj, and the ordering of the asymmetry in this region simply reflects 
the ordering of the different AG(?7) at large rj. Due to the relatively large asymmetry 
of up to 15%, such a measurement could provide vital information on the leading order 
polarized gluon distribution for 0.1 < x < 0.35 and ~ ^j/v-- 

In the foreseeable future, the HERA collider may be able to accelerate polarized pro- 



tons This would offer the unique opportunity of measuring the polarized structure 
function gi{x, Q^) far beyond the range of present fixed target experiments. In order 
to judge the quality of the new information gained from such a measurement, an order- 
of-magnitude estimate of the statistical errors is crucial [JT|. The structure function 
gi{x, Q"^) is extracted from a measurement of the asymmetry 

^ll(x,Q^) = AeA,D(y)||^^, (24) 
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where Ag (Ap) denote the polarizations of the electron (proton) beam and 



D 



(25) 



2il-y){l + Rp) + y^ 

is the kinematical depolarization of the photon. The statistical error on xgi{x^Q'^) is 



therefore ||41|| : 



1 



2(1-2/) 



;Fl{x,Q') 



l + (l-y)2- 

X I I [dV("'^P°i-V(da;dQ2)]dxdQ2') Jl - A||(x, Q^), (26) 



where the unpolarized differential cross section is integrated over the bin used in the 
measurement. 

To study the accuracy of a measurement of gi{x, Q"^) in the collider mode of HERA, 
we have evaluated the above expression using the unpolarized parton distribution set 
MRS (A') from Ref. with the next-to-leading order polarized parton distribution set 
A described in Section 2. We apply the following cuts to the HERA phase space ||4^ : 



0.1 < y < 0.95, Q > 1 GeV , Q^' < 176° and E^,' > 5 GeV, and consider two scenarios 
for the beam energies: 



[a) 
(b) 



300 GeV 
150 GeV 



with 



27 AA GeV, Ep 



820 GeV, 



with E, = 18.75 GeV, E„ = 300 GeV. 



(27) 



The expected errors for an integrated luminosity of 60 pb at Ag = Ap = 0.8 in these two 
scenarios are shown in Fig. 6. We take two bins per decade in x and only one bin in Q^. 



These results are consistent with the leading-order estimates of [^, bearing in mind the 
different cuts applied. In particular, we estimate smaller errors in the small-x region, as 
we apply no cuts on the hadroinc final state. 

The average values probed at HERA range from 2 GeV^ (x ^ 6-10-^) to 1060 GeV^ 
(x ^ 0.05) for 7i = 300 GeV and from 1.7 GeV^ (x ^ 2 ■ 10"^) to 270 GeV^ (x ^ 0.05) 
for ^/s = 150 GeV. For reference, we also show in Fig. 6 a parametrization of gfix, Q"^) 
obtained from the NLO set A. To illustrate the impact of a measurement at HERA, we 
include the three lowest datapoints reported by the SMC experiment |11|, corresponding 



to values around 1.5 GeV^. It is apparent that a measurement at lower beam energies 
will yield data of higher statistical quality. In contrast, the higher beam energies yield a 
measurement at smaller x. 

A measurement of gi{x, Q^) at the HERA collider will evidently not provide a large 
number of precision data on the x and dependences of this structure function. Hence, 
it will not provide sufficient information for an indirect determination of the polarized 
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gluon density. The important physics result of such a measurement is the determination 
of the small-x behaviour of gi{x, Q"^), for which various, significantly different predictions 
exist |2^. It is important to stress that a measurement of gfix, Q^) at small x will reduce 
the experimental uncertainty on the Ellis- Jaffe sum rule. The impact of the small x 
region can be easily seen in Fig. 6, as the Ellis- Jaffe sum rule is proportional to the area 
enclosed by xg^^x, Q^) and the x axis. 

A final remarkable point on the measurability of gi{x, Q"^) at HERA is the impact on 
the statistical error of the minimum cut on y. As the photon depolarizes for small y, a 
small cut on y (such as y > 0.01) diminishes the statistical accuracy, even though more 
data are included in the bin. We find that a minimal cut on y between 0.1 and 0.2 yields 
the optimal accuracy. 



4 Summary and Conclusions 

We have performed leading and next-to-leading order QCD fits to the world data on 
the gi polarized structure function measured with proton, neutron and deuteron targets. 
The quality of the fits is excellent, and from them we obtain sets of polarized parton 
distributions!^ which can be used for further phenomenological analyses. The experimental 
precision is highest for the proton and deuteron data, and together these constrain the 
shapes of the valence u and d distributions. The sea quark and gluon distributions are still 
largely undetermined. There is a weak constraint on the overall size of the former, but 
almost no information at all on the flavour decomposition of the sea. Following our earlier 
work, we have presented three qualitatively different gluon distributions, characterized by 
different behaviours at large x. 

It seems that further progress in the determination of polarized parton distributions 
will require measurements other than gi at fixed-target experiments. Inelastic J/ip pho- 
toproduction in lepton-hadron scattering has already proved capable of measuring the 
unpolarized gluon distribution, and a similar measurement with a polarized beam and 
target would provide the first real information on AG'(x, Q^). A measurement of the 
dependence of Q^) can help to estimate the overall amount of gluon polarization 

at higher and medium x. Althogh present experiments cannot provide such a measure- 
ment due to low statistics or too low values of Q^, first information can be expected 
from the next generation of fixed target experiments. We have also studied the accu- 
racy of a measurement of 5'i(a;,Q^) at the HERA collider with polarized beams. Such a 
measurement will help to discriminate between various predictions for the behaviour of 
gi{x, Q^) at small x, and will reduce the uncertainty on the EUis-Jaffe sum rule from the 
extrapolation to small x. 

^Thc FORTRAN code for the various sets described in this paper is available by electronic mail from 
T . K . GehrniannQdurham .ac.uk 
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Figure Captions 

[1] Leading order polarized parton distributions as described in the text at Qg = 4 GeV^ 
compared to the unpolarized distributions of [T^ . 

Next-to-leading order polarized parton distributions as described in the text at 
Qo = 4 GeV^ compared to the unpolarized distributions of 

Structure function measurements Q^) of the proton, deuteron and neutron 

compared to the next-to-leading order predictions obtained using Gluon A. 

dependence of Ai''^{x, Q'^) in next-to-leading order using the Gluons A, B and 
C. To illustrate the sensitivity of current experiments, we show the x = 0.035 
datapoints from the recent SMC and E143 measurements. 

Predictions for the asymmetry A{z) (Eq. (p2D ) in polarized J/ip photoproduction 
at = 45 GeV, using the different leading order polarized gluon distributions. 

Expected errors for a measurement of gi{x, Q"^) at the HERA collider with ^/s = 
300 GeV and y/s = 150 GeV. 
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Fig. 4 
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